Electronic Properties of Two-Dimensional Carbon 
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We present a theoretical description of the electronic properties of graphene in the presence 
of disorder, electron-electron interactions, and particle-hole symmetry breaking. We show that 
while particle-hole asymmetry, long-range Coulomb interactions, and extended defects lead to the 
phenomenon of self-doping, local defects determine the transport and spectroscopic properties. Our 
results explain recent experiments in graphitic devices and predict new electronic behavior. 
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Carbon with sp 2 bonding has many allotropic forms 
with different dimensionality: three-dimensional graphite 
0, one-dimensional nanotubes Q, zero-dimensional 
fullerenes 0, and two-dimensional Carbon, also known 
as graphene. Graphene can be considered the materia 
prima for the other forms of Carbon that can be obtained 
from it either by stacking (graphite), wrapping (nan- 
otubes), or creation of topological defects (fullerenes). 
Hence, the electronic response of many Carbon based 
systems depends fundamentally on the basic physics of 
graphene. Recent experiments in graphene-based devices 
have shown that it is possible to control their electrical 
properties, such as carrier type (particle or hole), by the 
application of external gate voltage 0, 0| . These experi- 
ments not only have opened doors to carbon-based nano- 
electronics but also pose new questions on the nature of 
the electronic properties in these systems. 

Being two-dimensional, true long range positional or- 
der of the Carbon atoms is not possible at any finite tem- 
perature, since thermal fluctuations introduce topologi- 
cal lattice defects such as dislocations (the Hohenberg- 
Mermim- Wagner theorem). Furthermore, because of the 
particular structure of the honeycomb lattice, the dynam- 
ics of lattice defects in graphene belongs to the generic 
class of kinetically constrained models 0, @ , where de- 
fects are never completely annealed since their number 
is only weakly dependent on the annealing time [|J. In- 
deed, defects are ubiquitous in carbon allotropes with 
sp 2 coordination, as confirmed by recent experiments ||. 
Furthermore, lattice defects induced by proton bombard- 
ment have been correlated to the appearance of mag- 
netism in graphitic samples 0, indicating the interplay 
between electron-electron interactions and disorder. 

Besides dislocations, graphene can sustain other types 
of extended defects such as disclinations, edges, and 
micro-cracks. It is known that certain edges (such as the 
zig-zag edge) lead to the appearance of localized states at 
the Fermi level Other defects such as pentagons 

and heptagons (lattice disclinations) also admit localized 
states [12j. Furthermore, tight-binding calculations show 
that a combination of a five-fold and seven-fold ring (a 
lattice dislocation) or a Stone- Wales defect (made up of 



two pentagons and two heptagons) also lead to a finite 
density of states at the Fermi level [la. Il4l ll5| . The pres- 
ence of states at the Fermi level generated by defects can 
be tracked down to the particular electronic structure of 
graphene. Each 7r-orbital contributes with one electron (a 
half-filled band) and the low-energy physics is described 
by a two-dimensional Dirac equation with a vanishing 
electronic density of states at the Fermi level. The van- 
ishing of the density of states has two very important 
consequences: the enhancement of the electron-electron 
interactions because of the absence of electronic screen- 
ing, and the formation of defect states at the Fermi level. 
In this respect, our work complements the extensive lit- 
erature on defects on electronic systems described by the 
two dimensional Dirac equation 16]. 

In the present work, we show that unscreened Coulomb 
interactions, particle- hole symmetry breaking and de- 
fects play a fundamental role in the electronic properties 
of graphene. In particular, graphene presents the phe- 
nomenon of self-doping, that is, charge transfer between 
extended defects and the bulk. In this case, depending 
on particle-hole symmetry breaking, graphene samples 
can either have electron or hole pockets. While self- 
doping derives from extended defects and controls the 
number and type of charge carriers, transport properties 
also depend on point-like defects such as vacancies, and 
surface ad-atoms. As we show in what follows, magneto- 
transport properties, such as Shubnikov-de Hass oscilla- 
tions and quantum Hall effect (QHE), as well as spec- 
troscopic quantities, such as quasiparticle lifetimes and 
infrared reflectivity, can be completely explained within 
this framework. Besides explaining published experimen- 
tal data, we also make new experimental predictions that 
can be used as tests of our theory. 

The bulk Hamiltonian of graphene can be written as: 
H = Ho + Hv + T~tu, where, 



wo = -t £ 4 



<r;«i,j» 



where Cj, ff (c\ a ) annihilates (creates) electrons at the 
site R; with spin a [a =T,I), t and t' are the near- 



est neighbor and next-nearest neighbor hopping ener- 
gies, respectively. For k — > the electronic dispersion 
is e k w 3i' + v F \k\ + (9i'a 2 |k| 2 )/4. Notice that t' does 
not change the structure of the extended electronic wave- 
functions, that remain described by the two-dimensional 
Dirac equation with a Fermi-Dirac velocity v-p — (3ta)/2 
(a is the lattice spacing). Tiy describes the disorder, 



(2) 



where Vj is the potential strength (for vacancies, Vj 
oo at the vacant site and is zero, otherwise), and rij 
^2i a Cj a Cj t o-. The electron interactions read, 



* ,3 



(3) 



e is the electric charge, and eo is the dielectric constant. 

Standard derivations of localized states near defects 
assume electron-hole symmetry, that is, t' — 0. In 
this case, the localized states lie exactly at the Fermi 
level. In order to ensure charge neutrality, these states 
are half filled, and there is no charge transfer between 
extended and localized states. Electron-hole symmetry 
is broken by terms in the Hamiltonian that allow for 
hopping between sites in the same sublattice, that is, 
when t' 0. In this case the band of localized states 
near an extended defect acquires a bandwidth of order 
t' and is shifted from Fermi energy. Using and 
we have studied in the Hartree approximation the in- 
duced charge transfer between the edge and the bulk 
in a long graphene ribbon of width W. In Fig.pQ we 
present the results for the induced electrostatic poten- 
tial, charge density, and charge transferred from a zig- 
zag edge to (delocalized) bulk states. We have studied 
zig-zag edges as long as O.lfi m in graphene ribbons with 
up to 10 5 atoms. As one can clearly see in Fig.pQ, self- 
doping effects are suppressed as the width of the edge 
increases. The Coulomb energy at the edge induced by 
a constant doping S per Carbon atom is approximately 
~ (Se 2 /a)(W/a). The charge transfer is arrested when 
the potential shifts the band of localized states to the 
Fermi energy, that is, when t' ps (e 2 /a)(W/a)6. Hence 
the self-doping is S ~ (f a 2 ) / (e 2 W) . The numerical re- 
sults are consistent with this estimate. Note that the sign 
of the induced charge depends on the sign of t'. The ex- 
tended defects lead to a scattering rate, T^ tic ~ vf/W. 
The pocket induced by the self-doping has a Fermi en- 
ergy of £f ~ (hvY^/S)/a ~ tyj (t 1 a 2 ) / (e 2 W) and thus, 
^"elastic ^ ep m wide ribbons. Hence, disorder due 
to extended defects does not smooth out the details of 
the small Fermi surface induced by the self-doping. For 
t = 2.7 eV Q, H = 0.2* ps 0.54 eV (see below), and 
a = 1.42A, we get 6 ~ 10~ 5 — 10~ 4 per Carbon atom in a 
system with domains of typical dimensions W ~ 0.1 — 1 fj, 
m, consistent with experiments |4|. 
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FIG. 1: Top: Electrostatic potential (dashed line) and elec- 
tronic charge (continuous line) as a function of position along 
a long graphene ribbon terminating on a zig-zag edge. En- 
ergies are given in units t, and distances in units of a. We 
have used t' = 0.2t, and e 2 /(eota) = 0.5. The inset shows the 
details of the electronic density near the edge. Bottom: Total 
additional charge per C atom in extended states as function 
of the ribbon width, W. 



We can use the ribbon geometry described above to 
analyze the properties of graphene in the integer QHE. 
A magnetic field B changes the phases of the hopping 
terms in the Hamiltonian leading to the appearance 
of bands of degenerate Landau levels. In the continuum 
limit the Landau levels have energy e n = -iiv-pl^ 1 y/n, 
where n is a positive integer, 1-q = yj &o/B is the cy- 
clotron radius, and <I>o = h/e is flux quanta. The Hall 
conductance can be obtained from the number of bulk 
states that cross the Fermi level due to the presence of 
an edge The calculated energy levels of a graphene 
ribbon in a magnetic field are shown in Fig. [2]. Notice 
that the momentum along the ribbon fixes the position of 
the levels. When the position approaches the edges, the 
positive energy levels rise in energy as for electron-like 
Landau levels in the 2D electron gas, while the negative 
energy levels behave in a hole-like fashion. There is a 
zero energy mode that splits into a set of electron and 
hole-like levels. Finally, the localized states at the edges 
are quite insensitive to the applied magnetic field. If we 
fix the chemical potential above (below) the zero mode 
energy, all bulk electron-like (hole-like) levels of lower 
energy give rise to crossings, and contribute to the Hall 
conductance. Each Landau level is doubly degenerated 
since the honeycomb lattice gives rise to two incquivalcnt 
sets of Dirac fermions . Hence, the Hall conductance 
arising from the bulk modes is: 



<7 Ha ii = 2(2N + l)e 2 /h , 



(4) 



and thus the Hall conductivity is quantized in odd num- 
ber of 2e 2 /h (the factor of 2 comes from the spin degen- 
eracy since cyclotron energy, huj c = y/2vFh/lB, is much 



3 



larger than the Zeeman energy, g^igB, where g sa 2 and 
[iB the Bohr magneton). This analysis neglects correc- 
tions due to the additional band of localized states at the 
edges induced by the structure of the ribbon. 





FIG. 2: Energy levels of a graphene ribbon of width 600 a 
with zig-zag edges as function of the momentum parallel to 
the edges. Energy in units of t, and t' — 0.2t. Momentum 
in units of l/(2\/3a). Top: zero magnetic field; Centre: a 
magnetic flux per C atom of $ = 0.00025 <E>o; Bottom: <3> = 
0.0005 $o- 
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FIG. 3: (a): Density of states as function of vacancy con- 
centration (D is a high energy cut-off); (b): As in (a), in the 
presence of an applied field. Dashed lines mark the positions 
of the Landau levels in the absence of disorder; (c) and (d): 
<t(/x) as function of the Fermi energy p. t' = in all cases. 

While extended defects induce self-doping they do not 
change significantly the electronic properties. However, 
point defects such as vacancies or ad-atoms can change 



the transport properties even if they do not lead to self- 
doping. It is known that weaker types of disorder do not 
change the semi-metallic nature of a system described 
by the Dirac equation 0|. Nevertheless, point defects 
in the unitary limit can lead to a sharp resonance at 
the Fermi energy ^(|. We describe the effect of vacan- 
cies in the transport properties using the Coherent Phase 
Approximation (CPA), that gives a good description of 
the spectral and transport properties of graphene, except 
perhaps within a small region near Fermi energy where 
electronic localization becomes important [2(|- Within 
CPA, the one-particle spectral function is written as: 



A(k, u) = Im [u - Scpa(w) - e k ] 



(5) 



where electron self-energy, Ecpa(w), has to be deter- 
mined self-consistently from 0]) and ©. The electronic 
density of states, p(u>) = ^k^C 4 ' 1 ^)' as function of ap- 
plied magnetic field and impurity concentration, in the 
continuum limit, is shown in Fig.[3j(a)-(b). Unlike for 
weaker forms of disorder, vacancies induce a finite density 
of states at the Fermi level. Furthermore, ImScpA(^) is 
finite at the Fermi level, and is monotonically increasing 
as u> — > indicating that vacancies have a strong effect 
on the Dirac fermions. This function can be considered 
an intrinsic linewidth as measured in ARPES. Besides 
the effect of disorder, the lifetime of quasiparticles has 
also an intrinsic contribution from interaction effects as- 
sociated with eq. © giving a contribution e 2 /vp |o;|[l9|. 
which is monotonically decreasing as lj — > 0. Therefore, 
the final result, as shown in Fig.^J], predicts that that 
the total linewidth of quasiparticles shows a minimum. 




FIG. 4: |S"(oj)|, as a function of energy to. Inset: A(k,ui). 
The frequency dependent conductivity can be written 



as: 



a(iu,T) = 2e 2 /(hir) / deK (e, u)[f(e, T)-f(e+u, T)]/u 

J — OO 

(6) 

where K(e, ui) is a conductivity kernel, and /(e, T) is the 
Fermi-Dirac distribution function. At low temperatures, 
the d.c. conductivity, c(/i) = <j(lo — 0,T), is approxi- 
mately proportional to K(fi,uj = 0), and can be experi- 
mentally measured by the application of a gate voltage, 
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V, to a graphene plane (in the presence of a voltage V 
the chemical potential shifts from // to fx + eV) . In Fig. 
[3] (c)-(d) we show a(/i). A non-conventional ^ depen- 
dence is found when compared to other 2D electronic 
systems, in agreement with experiments 0. As the tem- 
perature rises, the d.c. conductivity is found to increase 
with temperature, as the thermally excited carriers con- 
tribute (as in the case of a narrow gap semiconductor) 
and is also in agreement with the experimental data Q. 
When /i + eV — the low temperature conductance per 
plane has a universal value, independent of disorder and 
magnetic field |2l|, <7dc = 4e 2 / (71-/1). A similar universal 
behavior was predicted for d-wave superconductors \22\ . 
Because of particle-hole symmetry breaking, t' 7^ 0, we 
expect an asymmetry in the Shubnikov-de Haas oscilla- 
tions as function of gate voltage. We can quantify this 
asymmetry by looking at Landau levels j* and —j* — 1 
with j* ^> 1. It is possible to show Hj| that: 

\t'/t\ « V2/12(f + Z/2)-\f + l)-^ 2 (£ B /a) . (7) 

For j* w 6 for B = 12 T (£ B = 52 A) we find \t'/t\ w 
0.2, which is the value used in this work. 

In Fig.0 we show cf(uj) in the presence of impurities 
and finite magnetic field. The existence oscillations in 
the conductivity at low temperatures is clearly seen. An 
analysis of the transitions to the different peaks shows 
that those involving the j = Landau level are sup- 
pressed and that the energy of the lowest Landau levels 
is significantly shifted by the disorder. 
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FIG. 5: cr(w} for different magnetic fields and ni = 0.001, 
using the parameters given in the text. The lower right graph 
shows a scaling of the curves as function of iv/s/B. 

We have analyzed the influence of lattice defects on 
the electronic properties of graphene. Our results show 
that: (?) Extended defects lead to the existence of lo- 
calized states, or resonances, near the Fermi level of the 
pure system. In the absence of electron-hole symmetry, 
these states lead to self-doping, and to the existence of 



electron or hole Fermi pockets with 10~ 4 — 1CP 5 electrons 
per unit cell for domains of sizes 0.1 — 1/im. (ii) An inte- 
ger QHE with a quantized Hall conductivity in odd values 
of 2e 2 /h. (in) Point defects lead to a finite inverse scat- 
tering time in undoped graphene, which decreases with 
increasing energy. As the inelastic contribution rises lin- 
early with temperature, a minimum in the lifetime ob- 
served in ARPES experiments is predicted, (iv) The d.c. 
conductivity in undoped graphene has a universal value 
at low doping and temperatures. In undoped samples, it 
rises with temperature, as the number of thermally acti- 
vated carriers increases, (v) The optical conductivity in 
a magnetic field is modulated due to Landau levels. 
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